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1 | INTRODUCTION 


Abstract 

Increased evidence of porcine deltacoronavirus (PDCoV) causing diarrhoea in pigs 
has been reported in several countries worldwide. The virus has currently evolved 
into three separated groups including US, China and Southeast Asia (SEA) groups. In 
Vietnam, PDCoV was first reported in 2015. Based on phylogenetic analyses of spike, 
membrane and nucleocapsid genes, it is suggested that Vietnam PDCoV is chimeric 
virus. In the present study, we retrospectively investigated the presence of PDCoV in 
Vietnam and the full-length genomes of six PDCoV isolates identified in 2014-2016 
were further characterized. The results demonstrated that Vietnam PDCoV was first 
detected as early as 2014. All six Vietnam PDCoV are in the SEA group and further 
divided into two separated subgroups including SEA-1 and SEA-2. Vietnam PDCoV 
in SEA-2 was closely related to Thai and Lao PDCoV. Recombination analysis demon- 
strated that three isolates in SEA-1 were a chimeric virus of which P12_ 14 VN_0814, 
the first Vietnam isolate, and US PDCoV isolates were major and minor parents, re- 
spectively. The recombination was further evaluated by phylogenetic construction 
based on 3 recombinant fragments. The first and third fragments, closely related 
to P12 14 VN_0814, were associated with ORFia/1b and N genes, respectively. 
The second fragment, associated with S, E, and M genes, was closely related to US 
PDCoV isolates. High antigenic and hydrophobic variations were detected in $1 pro- 
tein. Three-day-old pigs challenged with the chimeric virus displayed clinical diseases 
and villus atrophy. In conclusion, Vietnam PDCovV is genetically diverse influenced by 
an external introduction from neighbouring countries. The chimeric Vietnam PDCoV 


can induce a disease similar to Thai PDCoV. 
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2010). PDCoV is the causative agent of a newly emerged enteric 


disease in pigs characterized by watery diarrhoea, dehydration and 


Porcine deltacoronavirus (PDCoV) is an enveloped, positive-sense, death with 40%-80% mortality following an outbreak (Jung, Hu, & 
single-stranded RNA virus belonged to the genus Deltacoronavirus, Saif, 2016). PDCoV infects enterocytes of the small intestine lead- 
family Coronaviridae, order Nidovirales (Woo, Huang, Lau, & Yuen, ing to villous atrophy and malabsorption (Hu et al., 2015; Jung et 
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al., 2015). Clinical signs caused by PDCoV resemble those with two 
viruses In the genus Alphacoronavirus, family Coronaviridae, includ- 
ing porcine epidemic diarrhoea virus (PEDV) and transmissible gas- 
troenteritis virus (TGEV) (Jung et al., 2015), but severity of clinical 
disease is lesser and piglet mortality is lower. PEDV and TGEV infect 
pigs of allages and cause 100% morbidity, with 80%-100% mortality 
rates in piglets. The higher severity of PEDV and TGEV infection in 
pigs indicates that the two viruses are more suitable for infecting 
pigs. On the other hand, PDCoV may not be perfectly adapted to 
pigs, and the virus may be evolving towards better ability to infect 
a specific host. Recently, an experimental study in calves suggested 
that calves are susceptible to infection by PDCoV but not PEDV 
(Jung, Hu, & Saif, 2017). 

The full-length genome of PDCoV is approximately 25 kb in 
length and contains seven essential genes, including open reading 
frame 1a/1b (ORF1a/1b), spike (S), envelop (E), membrane (M), and 
nucleocapsid (N), and specific accessory genes including non-struc- 
tural protein 6 (Nsp6) and Nsp7, flanked by a 5’- and 3’-untranslated 
region (UTR) (Woo et al., 2012). ORF1a and ORF1b, occupying two- 
thirds of the genome, encode 2 overlapping viral replicase/transcrip- 
tase polyproteins, 1a and 1b, which are cleaved into 15 Nsps, namely 
Nsp2-16 (Zhang & Yoo, 2016). The S gene encodes S protein which 
consists of two domains called $1 and S2 domain. The $1 domain 
plays an important role in binding to specific host cell receptors and 
contains a neutralizing epitope. The S2 domain functions in mem- 
brane fusion (Shang et al., 2018). The E and M genes are transmem- 
brane proteins associated with viral envelope formation and virus 
release (Woo et al., 2010). On the other hand, the functions of the 
N gene are associated with viral RNA replication and pathogenesis 
(Lee & Lee, 2015). 

PDCoV was first detected in Hong Kong (HKU15-44 and 
HKU15-155 isolates) in 2012 (Woo et al., 2012). Although it was first 
detected, there was no association with clinical disease. The first ev- 
idence of PDCoV causing a disease was first evident in Ohio, USA 
in 2014. Soon after the first emergence in Ohio, USA, the virus was 
then reported in 18 states (Jung et al., 2016). The rapid widespread 
of PDCoV was evident when the virus was subsequently reported 
in China, South Korea, Thailand, Laos, and Vietnam in 2015 (Le et 
al., 2018; Lee & Lee, 2014; Madapong et al., 2016; Marthaler, Jiang, 
Collins, & Rossow, 2014; Saeng-Chuto, Lorsirigool, et al., 2017; Song 
et al., 2015; Wang, Byrum, & Zhang, 2014a, 2014b). Although the 


Year Regions No. farms Total samples +ve PDCoV 
2016 Middle 2 6 0 
South 2 4 6: 
2015 South 14 2° 6 
2014 South 21 47 2 
203 South S Q 0 
2012 South 0 
2011 South 4 7 O 
Total 49 108 11 


+ve PEDV 


87 


first detection of PDCoV in those countries was in 2014-2015, the 
retrospective investigation of intestinal samples demonstrated that 
the presence of PDCoV in China, USA and Thailand was as early as 
2004, 2013 and 2013, respectively (Saeng-Chuto, Stott, et al., 2017; 
Sinha, Gauger, Zhang, Yoon, & Harmon, 2015). 

Presently, PDCoV has evolved into three separated groups, in- 
cluding US, China and Southeast Asia (SEA). PDCoV isolates from 
Thailand and Laos were clustered in a novel PDCoV group, SEA. The 
SEA group is genetically distinct from US and China PDCoV groups 
(Madapong et al., 2016; Saeng-Chuto, Lorsirigool, et al., 2017; Saeng- 
Chuto, Stott, et al., 2017). In Vietnam, PDCoV was first detected 
in 2015 (Saeng-Chuto, Lorsirigool, et al., 2017). Genetic analyses 
based on S, M, and N genes demonstrated that Vietnam PDCoV iso- 
lates were genetically different from the SEA group (Saeng-Chuto, 
Lorsirigool, et al., 2017). The phylogenetic analyses based on S and 
M genes demonstrated that Vietnam PDCovV is clustered in the SEA 
group. However, phylogenetic analysis based on the N gene demon- 
strated contradictory results, since Vietnam PDCoV was clustered 
in the US group. The findings implied that Vietnam PDCoV could 
potentially be a recombinant virus. Whether or not this genotype 
of PDCoV isolates predominantly existing in Vietnam remains not 
known. We therefore conducted a retrospective study investigating 
the presence of PDCoV in intestinal samples collected from pigs less 
than 5 days old experiencing severe diarrhoea, in Vietnam from 2011 
to 2016. The full-length genome sequences of the detected Vietnam 
PDCoV were further characterized on the basis of heterogenicity, 
and recombination analyses together with virus infectivity in cell cul- 


ture and pathogenicity in piglets were performed. 


2 | MATERIALS AND METHODS 


2.1 | Samples and viruses 


One hundred and eight intestinal samples collected from 49 swine- 
herds with clinical diarrhoea outbreaks in different regions in Vietnam 
(Table 1) were collected for viral pathogens causing diarrhoea during 
2011-2016. All intestinal samples were individually minced using an 
Axygen™ Tissue Grinder (Thermo Fisher Scientific) and resuspended 
in 1X phosphate-buffered saline (PBS; 0.1 M, pH 7.2). The suspen- 
sions were centrifuged at 10,000 x g for 10 min at 4°C, and the 


supernatant filtered through a 0.45-um nylon membrane (Corning 


TABLE 1 _ Intestinal sample number 
that was collected from different regions 
in Vietnam in 2011-2016 


+ve TGEV 
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Inc.). Virus was isolated, and plaque purification was performed as 
described in Data S1. 


2.2 | Virus detection and full-length 
genome sequencing 


Viral RNA was extracted from the supernatant using Nucleospin 
Viral RNA Extraction Kit (Macherey-Nagel Inc.) and converted to 
cDNA using M-MuLV Reverse Transcriptase (New England Biolabs 
Inc.). All samples were tested with the specific primer for S gene of 
PEDV and specific primer for N gene of TGEV as in previous studies 
(Lee & Lee, 2014; Park et al., 2007; Song et al., 2015; Temeeyasen et 
al., 2014; Wang, Byrum, & Zhang, 2014a; Woo et al., 2010). PDCoV 
was detected by PCR using specific primer for M and N genes as pre- 
viously described (Wang, Byrum, & Zhang, 2014b). The PCR prod- 
ucts were electrophoresed at 100 V for 30 min on a 1% agarose gel, 
and the gel was stained with RedSafe™ nucleic acid staining solution 
(INtRON Biotechnology Inc.) and examined under a UV light. The 
PCR products were purified using Nucleospin Gel and PCR Clean-up 
kit (Macherey-Nagel Inc.), and sequencing was performed by First 
BASE Laboratories Inc. PDCoV-positive samples were subjected to 
whole-genome sequencing using 26 primers (Table $1). All specific 
PCR products were purified using Nucleospin Gel and PCR Clean-up 
kit (Macherey-Nagel Inc.), and sequencing was done by First BASE 
Laboratories Inc., using an ABI Prism 3730XL DNA sequencer. 


2.3 | Phylogenetic and genetic analyses 


To determine the genetic relationship between the six Vietnam 
PDCoV isolates and other PDCoV isolates, the nucleotide (nt) and 
amino acid (aa) sequences were aligned using the CLUSTALW pro- 
gram (Thompson, Higgins, & Gibson, 1994), together with 84 other 
PDCovV isolates available in GenBank (Table S2). The nt and aa iden- 
tities were determined using the sequence identity matrix function 
implemented in BioEdit software (Hall, 1999). The phylogenetic tree 
based on the full-length genome, open reading frame (ORF) 1a/1b, 
SEMN and S, M and N genes was separately built using the maximum 
likelihood (ML) method with 1,000 bootstrap replicates, and substi- 
tution models were selected by the best-fit substitution model func- 
tion implemented in MEGA software version 6.0 (Tamura, Stecher, 
Peterson, Filipski, & Kumar, 2013). TN93 + G+l model was used in 
phylogenetic trees based on full-length genome and SEMN gene 
construction. TN93 + G model was used in phylogenetic trees based 
on ORF1a/1b and S genes construction. K2 + G and K2 + G+l models 
were used in the construction of phylogenetic trees based on M and 


N genes, respectively. 


2.4 | Recombination analysis and phylogenetic 
reconstruction based on recombination fragments 


Recombination of six Vietnam PDCoV isolates was further analysed 
together with the 84 reference PDCoV isolates using automated 
RDP, GENECONV, BOOTSCAN, MaxChi, CHIMAERA and SISCAN 


methods implemented in RDP4 software (Martin, Murrell, Golden, 
Khoosal, & Muhire, 2015). Major and minor parents and breakpoints 
also were detected, and recombination fragments were generated 
with the RDP4 software (Martin et al., 2015). 

Phylogenetic trees based on recombinant fragments were sep- 
arately constructed using the ML method with 1,000 bootstrap 
replicates, and substitution models were designated by the best-fit 
substitution model function implemented in MEGA software version 
6.0 (Tamura et al., 2013). TN93 + G+l model was used in phyloge- 


netic tree based on all recombination fragments construction. 


2.5 | Sliding window, antigenicity and 
hydrophobicity analyses 


The sliding window, hydrophobicity and antigenicity were evaluated 
among the chimeric Vietnam PDCoV isolates in this study, and sus- 
pected major and minor parents of the chimeric Vietnam PDCoV, 
respectively (Figure 2). 

The nt sequence variation sites were predicted using dnaSP 
version 6 with 100 bp and a step size of 25 bp (Rozas et al., 2017). 
Antigenicity and hydrophobicity were also analysed using IEDB 
Analysis Resource (http://tools.immuneepitope.org/main) with the 
Kolaskar and Tongaonkar method and ExPASy-ProtScale (http://web. 
expasy.org/protscale) with the Kyte-Doolittle method, respectively. 


2.6 | Invitro infection of recombinant 
Vietnam PDCoV 


Viral growth curves of the three chimeric Vietnam PDCoV isolates 
and NT1_1215 isolate were compared. Each isolate was performed 
in quadruplicate. The NT1_1215 isolate is a Thai PDCoV isolated 
from a swineherd in Thailand experiencing a diarrhoea outbreak 
with piglet mortality at approximately 40%. The full-length genome 
is deposited in GenBank with accession number KX361345. 

The PDCoV-positive samples were subjected to virus isola- 
tion using swine testicular (ST) cells (ATCC CRL-1746). The viral 
isolation was approved by Institutional Biosafety Committee of 
Chulalongkorn University (CU-IBC) (IBC 1831062). In brief, ST cells 
were seeded in 125 flasks (Corning, NY, USA) containing Advanced 
Minimum Essential Medium (adv. MEM) supplemented with 10% 
heat-inactivated foetal bovine serum (FBS), 1% HEPES, 1% antibi- 
otic-antimycotic, 1% L-glutamine (Gibco) and 0.0001% pancreatin 
from porcine pancreas (Sigma-Aldrich) for each viral inoculation. 
After the cells reached 80% confluency, the cells were inoculated 
with each virus at a titre of 10° TCID,,/ml. The infected cells were 
incubated at 37°C and 5% CO... The viruses were harvested at dif- 
ferent time points (0, 12, 24, 36, 48, 60 and 72 hr post-inoculation 
(hpi.)) and stored at -80°C. 

Viral titre at each time point was determined. ST cells were cul- 
tured in 96-well plates (Corning, NY, USA) containing supplemented 
adv. MEM. Cells were inoculated with each virus at each time point 
at 10-fold serial dilutions (10°-10’). Cytopathic effect was deter- 


mined daily, and the growth curves were plotted. 
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2.7 | In vivo infection of recombinant 
Vietnam PDCoV 


Eighteen 3-day-old healthy piglets were randomly divided into three 
groups: (i) mock control group (n = 6), (ii) chimeric Vietnam PDCoV 
(P29_15_ VN_1215) isolate-infected group (n = 6), and (iii) NT1_1215 
isolate-infected group (n = 6). Each piglet was inoculated with 5 ml 
of each virus at a titre of 10° TCID,,/ml. Clinical signs were observed 
daily. Three piglets per group were euthanized at 3 and 5 days post- 
inoculation (dpi.), and five segments of the small intestine, including 
duodenum, proximal jejunum, middle jejunum, distal jejunum and 
ileum, were collected and fixed in 10% formalin. Haematoxylin and 
eosin (H&E) staining was performed, and villus/crypt ratio was cal- 
culated. Comparison between each group was performed with one- 
way analysis of variance, followed by Tukey's multiple comparison 
test implemented in GraphPad Prism 7 (GraphPad Software Inc.). 


p < .05 was considered significant. 


3 | RESULTS 


3.1 | PDCoV detection 


One hundred and eight intestinal samples collected from pigs with 
diarrhoea were assayed for the presence of three viral pathogens 
including PEDV, TGEV and PDCoV using PCR assays (Table 1). Of 
108 intestinal samples tested, none were positive for TGEV. Eleven 
(10.19%) and 87 (80.56%) samples were positive for PDCoV and 
PEDV, respectively. All PDCoV-positive samples were also positive 
for PEDV. PDCoV was detected in Vietnam as early as 2014, but the 
most detection was mainly in 2015-2016. Of 11 samples, only six 
samples including P12_14 VN_0814, P29_15 VN_1215, P30_15_ 
VN_1215, P1_16_VN_0116, P19_16_VN_0416, and P20_16_ 
VN_0416 isolates were able to characterize their full-length genome 


sequences and further genetic analyses. 


3.2 | Full-length genome characterization 


The full-length genome sequence of the six Vietnam PDCoV iso- 
lates, recognized as P12_14 VN_0814, P29_15 VN_1215, P30_15_ 
VN_1215, P1_16_ VN_0116, P19_16_VN_0416 and P20_16_VN_0416, 
were deposited in the GenBank database under accession numbers 
MH700628, KX998969, MH118333, MH118331, MH118332 and 
MH700629, respectively. The six isolates were varied in size, varying 
from 25,405-25,409 nt. The major difference is the deletion in nu- 
cleotides in ORF1a/1b region. The full-length genome sequence of 
P12 14 VN_0814 isolate was 25,406 nt in length excluding poly(A) 
tail. P29_15_VN_1215, P30_15_VN_1215 and P1_16_VN_0116 iso- 
lates were 3 nt longer in which their full-length genome sequences 
were 25,409 nt in length excluding poly(A) tail. The full-length genome 
sequence of P19_16 VN _ 0416 and P20 16 VN_0416 isolates were 
25,405 nt in length excluding poly(A) tail. All isolates were character- 
ized by the same gene order of the 5'-UTR, ORF1a/1b, S, E, M, Nsp6, 
N, Nsp7 and 3’-UTR. The length of each region is presented in Table S3. 


Pairwise nt and aa sequence identity values of each gene be- 
tween the six Vietnam PDCoV isolates together with two previously 
reported Vietnam PDCoV isolates, HaNoi6/2015 and Binh21/2015 
isolates, and other reference PDCoV isolates are presented in 
Table S4. Based on the full-length genome analysis, the P12_14_ 
VN_0814, P29_15 VN_1215, P30_15_VN_1215 and P1_16_ 
VN_0116 isolates were closely related to previously reported 
HaNoi6/2015 and Binh21/2015 isolates (98.7%-100% nt) than 
Thai and Lao PDCoV isolates (98.1%-98.4% nt), US PDCovV isolates 
(97.5%-98.1% nt), Japan PDCoV isolates (97.2%-98.1% nt), South 
Korean PDCovV isolates (97.2%-98.0% nt) and China PDCoV iso- 
lates (97.1%-98.3% nt). In contrast, P19_16 VN_0416 and P20 _16_ 
VN_0416 isolates were more closely related to the Thai and Lao 
PDCoV isolates (99.1%-99.8% nt identity) than the HaNoi6/2015 
and Binh21/2015 isolates (98.3%-98.4% nt), US PDCoV isolate 
(97.2%-97.5% nt), China PDCoV isolates (96.9%-97.7% nt), Japan 
PDCoV isolates (96.9%-97.3% nt) and South Korean PDCoV iso- 
lates (96.9%-97.3% nt). 

The ORF1a/1b gene of all six Vietnam PDCoV isolates owns aa 
deletion and insertion in their genomes compared to other PDCoV 
isolates. The deletions of two (4°°LK*°) and three (7°8PVG”®) aa 
were observed at positions 401 to 402 and 758 to 760, respectively. 
Interestingly, these 5 aa deletions in the six Vietnam PDCoV isolates 
were similar to other PDCoV isolates from SEA countries. The inser- 
tion of one aa (°*N) at position 51 of the S gene was found in all six 
Vietnam PDCoV isolates similar to US PDCoV isolates. Based on the 
E, M and N genes, all six Vietnam PDCoV isolates did not show any 


aa deletion or insertion compared to other PDCoV isolates. 


3.3 | Phylogenetic analysis 


Phylogenetic trees constructed based on the full-length genome 
and ORF1a/1b demonstrated that the six Vietnam PDCoV isolates 
were grouped in the SEA group together with Thai, Lao, and two 
previously reported Vietnam PDCoV isolates based on full-length 
genome and ORF1a/1b gene sequences. The SEA group was further 
evolved into two subgroups including subgroup SEA-1 and SEA-2 
(Figure 1a,b). P19_16_VN_0416 and P20 16 VN_0416 together 
with Thai and Lao PDCoV isolates were clustered in the subgroup 
SEA-1. In contrast, P12_14_ VN_0814, P29_15_VN_1215, P30_15_ 
VN_1215 and P1_16_VN_0116 were clustered in subgroup SEA-2 
(Figure 1a,b). Interestingly, PDCoV isolates in the subgroup SEA-2 
are from Vietnam only. 

Phylogenetic trees constructed based on SEMN and S genes 
demonstrated that P19_16 VN 0416 and P20 16 VN_0416, and 
P12 14 VN_0814 were clustered together with the SEA group, sub- 
group SEA-1 and SEA-2, respectively (Figure 1c,d). However, three 
of the six Vietnam PDCoV isolates, including P29_15_VN_1215, 
P30_15_VN_1215 and P1_16 VN_0116, were clustered together 
with the US group (Figure 1c,d). 

The phylogenetic tree based on N gene demonstrated that 
P12 14 VN_0814, P29_15_VN_1215, P30_15_VN_1215 and P1_16_ 
VN_0116 were clustered together with the China group, differing 
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FIGURE 1 Phylogenetic trees based on nucleotide sequences of full-length genome (a), ORF1a/1b (b), SEMN (c), S (d), M (e) and N (f) 
genes. The phylogenetic trees were constructed using the maximum likelihood (ML) method with 1,000 bootstrap replicates implemented 
in MEGA software version 6.0. Red dots represent the PDCoV isolates from Vietnam in this study. All phylogenetic trees, except M gene, 
consisted of 3 groups namely US, China, and Southeast Asia (SEA) groups. SEA group was divided into 2 subgroups namely SEA-1 and SEA-2 
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Tract of sequence with a recombinant origin 


~~ P12 14 VN_0814 (Major Parent) - USA/Minnesota159/2014 (Minor Parent) 
~~ P12 14 VN_ 0814 (Major Parent)— P29 15 VN_ 1215,P30 15 VN 1215,andP1 16 VN 0116 (Recombinant) 
— USA/Minnesota159/2014 (Minor Parent) - P29 15 VN_ 1215, P30_15 VN_1215,and P1_16 VN_0116 (Recombinant) 


FIGURE 2 Recombination plot between P12 _ VN_14 0814, USA/Minnesota159/2014 and three Vietnam PDCoV isolates: P29_15_ 
VN_1215, P30_15_VN_1215 and P1_16_VN_0116 isolates as the major parent, minor parent and recombinant, respectively 


from P19_16 VN_0416 and P20 16 VN_0416, which were clus- 
tered in the SEA group (Figure 1f). 


3.4 | Recombination analysis and phylogenetic 
reconstruction based on recombination fragments 


Three recombination fragments were generated. The first and 
second fragments are at nt 1-19,453 and 23,511-25,426 (second 
fragment) in which are locations of ORF1ia/1b and N genes, re- 
spectively. The third fragment is at nt 19,454-23,512 in which is 
located at S, E and M gene region (Figure 2). Recombination analy- 
sis suggested that three Vietnam PDCoV isolates including P29_15_ 
VN_1215, P30_15 VN_1215 and P1_16 VN_0116 isolates were 
recombination of the P12_14 VN_0814 isolate (major parent) at 


the first and third fragments, and USA/Minnesota159/2014 isolate 
(minor parent) at the second fragment (Figure 2). The phylogenetic 
trees based on the first and third fragments demonstrated that 
these three chimeric Vietnam PDCoV isolates were closely related 
to the P12 14 VN_0814, Hanoi/2015 and Binh21/2015 isolates 
(Figure 3a,c). However, the phylogenetic trees based on the second 
fragment suggested that that the three chimeric Vietnam PDCoV 
isolates were closely related to US PDCoV isolates (Figure 3b). 


3.5 | Sliding window, hydrophobicity and 
antigenicity analyses 


The full-length genome sequences of the five isolates includ- 
ing P29.15 VN_1215, P30_15.VN_1215, P1_16_VN_0116, 


FIGURE 3 _ Reconstruction of phylogenetic trees based on the nucleotide sequence of three recombinant fragments, including first 
(a), second (b) and third (c) fragments. The phylogenetic trees were constructed using the maximum likelihood (ML) method with 1,000 
bootstrap replicates implemented in MEGA software version 6.0. Red dots represent the PDCoV isolates from Vietnam in this study. Blue 
triangle represents P12_VN_14 0814 isolate (major parent). Green triangle represents USA/Minnesota159/2014 isolate (minor parent) 
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FIGURE 3. (Continued) 


P12 14 VN_0814 (major parent) and USA/Minnesota159/2014 
(minor parent) were used in the sliding window analysis. Five sites, 
designated R1, R2, R38, R4 and R5, demonstrated high degree of 
genetic variation between all five isolates. Three sites including R1 
(nt position 1,326-1,425), R2 (nt position 2,101-2,225) and R3 (nt 
position 5,601-5,775) were located within ORF1a/1b regions. Two 
sites including R4 (nt position 20,101-20,200) and R5 (nt position 
25,001-25,100) were located within the S region, and between the 
N and 3’-UTR regions (Figure 4). 

The aa of the S gene encoding the structural protein associated 
with binding and neutralizing epitopes were subjected for further 
antigenicity and hydrophobicity analyses. Six segments at positions 
aa 7-15, 36-44, 116-124, 134-139, 346-352 and 546-554 demon- 
strated a difference in hydrophobicity and antigenicity indices be- 


tween the five isolates (Figure 5). 


3.6 | Invitro infection with recombinant 
Vietnam PDCoV 


Viral growth curves of the chimeric Vietnam PDCoV (P29_15_ 
VN_1215) and Thai (NT1_1215) isolates in ST cells were compared. 
The two isolates had a similar pattern of viral growth curves, but 
virus titres were different (Figure 6). Cytopathic effect, character- 
ized by cell rounding and clumping, with both PDCoV isolates were 
first observed at 12 hpi with a virus titre of 40° TCID.,/ml (mean 
log,, TCID,,/ml + SD; 1 + 0). At 24 hpi, the titre of P29_15_VN_1215 
remained at 10! TCID,,/ml (mean log,, TCID.,/ml + SD; 1 + 0), while 
the titre of NT1_1215 increased to 10 TCID.,,/ml (mean log,, TCID./ 
ml + SD; 2 + 0). The titres of both PDCoV isolates slowly rose from 
24 to 48 hpi. At 36 and 48 hpi, the titres of P29_15 VN_1215 were 
107 (mean log,, TCID.,/ml + SD; 2 + 0.25) and 10* TCID,,/ml (mean 
log,, TCID,,./ml + SD; 4 + 0.25), respectively, which were lower than 
those of NT1_1215 at 10° (mean logj5 TCID;4/ml SD; 3 £0.25) and 
10° TCID,,/ml (mean log,, TCID,,/ml + SD; 5 + 0.25), respectively. 
The titres of both PDCovV isolates rapidly declined to 10* TCID,,/ 
ml at 60 hpi (mean log,, TCID,,./ml + SD; 1 + 0.25). However, at 72 
hpi., the titre of NT1_1215 increased to 110+ TCID,,/ml (mean lOog,o 
TCID.,./ml + SD; 2 + O), but the titre of P29_15_VN_1215 remained 
at 10° TCID.,/ml (mean log,, TCID,,/ml + SD; 1 + 0). 


3.7 | In vivo infection with recombinant 
Vietnam PDCoV 


Pigs in the mock control group of which orally inoculated with culture 
media displayed no clinical disease throughout the study. In contrast, 
pigs orally inoculated with either Vietnam PDCoV isolates (P29_15_ 
VN_1215) or Thai PDCoV isolate (NT1_1215 isolate) displayed simi- 
lar level of clinical diseases associated with PDCoV including severe 
diarrhoea, vomiting, dehydration, weakness and lethargy. Clinical 
diseases were observed in both Vietnam and Thai PDCoV infected 
groups at 1 dpi. However, no piglets died before euthanasia. 

At 3 and 5 dpi, 3 piglets in each group were necropsied 


and five parts of small intestine including duodenum, proximal 
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FIGURE 4 _ Differences in full-length genome sequence between five isolates, including P29_15 VN_1215, P30_15_VN_1215, P1_16_ 
VN_0116, P12_VN_14_ 0814 (major parent) and USA/Minnesota159/2014 (minor parent), were analysed by sliding window method. The red 
arrows show high variation regions 
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FIGURE 5_ Hydrophobicity and antigenicity differences were found in S protein between the five isolates studied: P29_15 VN_1215, 
P30_15 VN_1215, P1_16_VN_0116, P12_VN_14 0814 (major parent) and USA/Minnesota159/2014 (minor parent) 


jejunum, middle jejunum, distal jejunum and ileum were col- piglets infected with P29_15 VN_1215 and NT1_1215 showed 
lected for further histopathological examination. H&E staining villous atrophy compared to that of the mock control group (G1) 


was performed. All H&E-stained small intestine segments of (Figure 7). 
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FIGURE 6 Viral growth curves of the Thai PDCoV (NT1_1215) (green line) and the chimeric Vietnam PDCoV (P29_15_ VN_1215) (red line) 
isolates. Statistical significance between groups annotated by same lower-case letters. p-value < 0.05 


Villus/crypt (V/C) ratios of piglets in each group were determined 
(Figure 8). At 3 dpi, the V/C ratio in the three jejunum segments of 
piglets infected with both PDCoV isolates was significantly lower 
compared to that of the mock control group (p < .05) (Figure 8b- 
d). At 5 dpi, the V/C ratio in middle and distal jejunum segments of 
piglets infected with both PDCoV isolates was significantly lower 
than that in piglets in the mock control group (p < .05) (Figure 8c,d). 
Although no significant difference was observed, the V/C ratio in all 
parts part in both infected groups trend to reduced when compared 


between 3 and 5 dpi within the groups (Figure 8a,c). 


4 | DISCUSSION 


In Vietnam, the first detection of PDCoV was reported in 2015 in the 
southern region of Vietnam (Saeng-Chuto, Lorsirigool, et al., 2017). 
Genetic analyses based on S, M and N genes discovered contradic- 
tory results suggesting that Vietnam PDCoV could be a recombinant 
virus. However, the recombination of the virus has not been con- 
firmed and the full-length genome was not attempt in the previous 
study. We therefore conducted a retrospective study investigating 
the presence of PDCoV in Vietnam, genetic diversity and character- 
ize the full-length genome of detected PDCoV. The isolated chimeric 
virus was further evaluated for its pathogenicity. 

In the present study, PDCoV was first detected in southern re- 
gion of Vietnam as early as 2014. The full-length genome sequences 
of six Vietnam PDCoV isolates were characterized. Phylogenetic 
analysis based on the full-length genome demonstrated that 
Vietnam PDCoV was genetically diverse. The six Vietnam PDCoV 
isolates were grouped in the SEA cluster together with the Vietnam 
(Hanoi/2015 and Binh21/2015), Thai and Lao PDCoV isolates which 
all have an amino acid deletion at positions 401 to 402 and 758 to 
760 located in the ORF1a gene. However, the six Vietnam PDCoV 


isolates were separated into two sub-clusters recognized as SEA-1 


and SEA-2. Two of six isolates were clustered in SEA-1 together with 
Thai and Lao PDCoV isolates while other four isolates were clus- 
tered in SEA-2 together with Hanoi/2015 and Binh21/2015 isolates. 

Due to a separated origin, the findings from the present study 
suggest there were at least 2 possible external introductions 
of PDCoV into Vietnam. As we mentioned earlier, the Vietnam 
PDCoV isolated in 2016 (P19_16_VN_0416 and P20_16_ VN_0416) 
shares high degree of genetic similarities and was clustered to- 
gether with Thai and Lao PDCoV in subgroup SEA-1, the virus 
isolates might be introduced from these neighbouring countries. 
Secondly, the external introduction of US-related PDCoV was 
also observed. This introduction seemed to influence the evo- 
lution of Vietnam PDCoV as evidenced by chimeric PDCoV in 
subgroup SEA-2. The Vietnam PDCovV isolates detected in 2014 
(P12 VN_14 0814) were clustered together with other Vietnam 
PDCoV isolates, HaNoi6/2015 and Binh21/2015, in subgroup 
SEA-2. The HaNoi6/2014 isolate was detected in Ha Noi province 
in the northern region of Vietnam since October 2015 and consid- 
ered as recombinant virus between Thai and China isolates (Le et 
al., 2018). This study shows the P12_VN_14 0814 isolate shares 
100% similarities with the HaNoi6/2014 isolate which means the 
P12 VN_14 0814 isolate could be also recombinant virus be- 
tween Thai and China isolates. Interestingly, three of other four 
isolates in subgroup SEA-2 were different. The P29_15 VN_1215, 
P30_15 VN_1215 and P1_16_VN_0116 isolates were confirmed 
to be chimeric virus which were a recombination between 
Vietnam PDCoV (P12_VN_14 0814), the major parent, and US 
PDCoV (USA/Minnesota159/2014), the minor parent. As previous 
reports, the recombination in RNA virus occurs when a host cell 
is infected with at least two viral genomes at the same time and 
affects host cell tropism, immune evasion, virulence and patho- 
genesis (Martin et al., 2011; Simon-Loriere & Holmes, 2011). The 
findings in this study suggested that there are at least four differ- 


ent genome sequences of PDCoV currently circulating in swine 
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FIGURE 7 Haematoxylin and eosin (H&E) staining of five segments of small intestine, including duodenum (a-c), proximal jejunum (d-f), 
middle jejunum (g-i), distal jejunum (j-l) and ileum (m-o). Villous atrophy was observed in all small intestine segments of piglets infected with 


the chimeric Vietnam PDCoV (P29_15_ VN_1215) (b, e, h, k and n) or the Thai PDCoV (NT1_1215) (c, f, i, land 0) isolates, compared to mock 
control group (a, d, g, j and m). Statistical significance between groups annotated by asterisk. p-value < 0.05 
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FIGURE 8 _ Villus to crypt (V/C) ratio of five segments of small intestine, including duodenum (a), proximal jejunum (b), middle jejunum (c), 


distal jejunum (d) and ileum (e) 


farms in Vietnam. The virus may be trying to evolve, adapt to pigs, 
emerge and cause a disease. It is remained to be seen whether 
each of four different genome sequence groups would influence 
the genetic evolution of each other and become a dominant group, 
or continuously recombine and evolve in the population, resulting 
in the emergence of a novel viral isolate. 

To analyse the nt variation between chimeric Vietnam PDCoV 
isolates and their parents, the sliding window method was used. The 
results showed that high nt variation was detected in five sites of the 
full-length genomes. Most nt variations occurred in the ORF1a/1b 
gene. One of the five sites was found in the S gene which encodes 
the S protein. The hydrophobicity and antigenicity variations of the S 
protein were then analysed. Hydrophobicity is associated with pro- 
tein-protein interaction (e.g. ligand of virus-host receptor binding) 
and can be changed by replacement with hydrophobic amino acids 
(Bosch, Zee, Haan, & Rottier, 2003; Tekewe, Connors, Middelberg, 
& Lua, 2016), and change in antigenicity can affect immune re- 
sponse (Group et al., 2012; He, Li, Heck, Lustigman, & Jiang, 2006). 
Differences in hydrophobicity and antigenicity were detected in six 
parts in the $1 domain (aa 1-552) among the five isolates. The S1 
domain plays an important role in host receptor binding and contains 
a neutralizing epitope (Shang et al., 2018; Woo et al., 2010). These 
results might support the proposition that the chimeric Vietnam 
PDCoV isolates are evolving towards improved host infection and 
immune response evasion. 

To investigate the infectivity of the chimeric Vietnam PDCoV 
(P29_15_VN_1215) isolate compared to the Thai isolate (NT1_1215), 


both virus isolates were propagated in ST cells. The viral growth 


curves were developed and compared to each other. The results 
showed that the growth curve of chimeric Vietnam PDCoV isolate 
(P29_15_VN_1215) in ST cells was relatively lower compared to the 
Thai isolate at 12-60 hpi. This finding suggested that the chimeric 
Vietnam PDCoV isolate (P29_15_VN_1215) had less infectivity in ST 
cells than the Thai isolate. 

To demonstrated whether or not this chimeric Vietnam PDCoV 
(P29 15 VN_1215) isolate can induce clinical disease, the virus 
was orally inoculated in piglets. Clinical diseases were compared 
to piglets that were orally inoculated with the Thai isolate. The re- 
sults demonstrated that both PDCoV isolates can induce clinical 
diseases and pathogenicity in piglets. V/C ratio of piglets infected 
with chimeric Vietnam PDCoV was reduced but not significant dif- 
ference compared to those of piglets infected with Thai PDCoV. 
The PDCoV infectivity results are in line with other studies. In 
2015, two PDCoV isolates, OH-FD22 and OH-FD100, were also 
shown to cause disease. The infected gnotobiotic pigs showed 
severe diarrhoea, vomiting and severe atrophic enteritis (Jung et 
al., 2015). Moreover, in 2018, piglets were challenged with PDCoV 
KNU16-07 isolate and showed similar clinical signs, namely diar- 
rhoea, vomiting and viral enteritis (Jang et al., 2018). The findings 
in our study confirmed that the chimeric Vietnam PDCoV infection 
could induce small intestinal enteritis. 

In conclusion, Vietnam PDCoV is genetically diverse influencing 
from the external introduction. Vietnam PDCoV is further evolved 
into 2 separated subgroups including SEA-1 and SEA-2. Vietnam 
PDCoV in SEA-2 was closely related to Thai and Lao PDCoV. 


Recombination analysis demonstrated that isolates in SEA-1 were a 
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chimeric virus. P12_14 VN_0814, the first Vietnam isolate, and US 
PDCoV isolates were major and minor parents, respectively. The chi- 


meric Vietnam PDCoV could induce clinical diseases. 
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